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We present an Anderson-type model Hamiltonian with exchange coupling between the localized 
spins and the confined holes in the quantum dots to study the ferromagnetism in diluted magnetic 
semiconductor (DMS) quantum dot arrays embedded in semiconductors. The hybridization between 
the quantum-confined holes in the DMS quantum dots and the itinerant holes in the semiconductor 
valence band makes hole transfer between quantum dots, which can induce the long range ferromag- 
netic order of the localized spins. In addition, it makes the carrier spins both in the DMS quantum 
dots and in the semiconductors polarized. The spontaneous magnetization of the localized spins 
and the spin polarization of the holes are calculated using both the Weiss mean field approximation 
and the self-consistent spin wave approximation, which are developed for the present model. 
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I. INTRODUCTION 

Diluted magnetic semiconductors (DMS) and a vari- 
ety of quantum nanostructures based on DMS materials 
have recently attracted much interest due to the novel 
physics and the potential application to the emergent 
field of spintronicsa The discovery of the ferromagnetism 
in the III-V based DMS materials Gai-ajMn^As and 
Ini-^Mn^As has made it possible to combine the mag=. 
netic and semiconducting properties in one material.H 
At high concentration of randomly distributed Mn 2+ 
ions doped in GaAs samples with high hole density, 
Gai^^Mn^As compounds exhibit ferromagnetism with 
transition temperature as high as 110 K at certain value 
of x. The ferromagnetism in the localized spins of impu- 
rity Mn ions is mediated by the itinerant holes through 
the p-d exchange coupling Vfi^WPpn ^ ne va lance-band 
holes and the localized spinsBlaQEl Bhatt et al. argues 
that for low carrier density and strong disorder limit, the 
ferromagnetism is mediated by the carriers in impurity 
bands.d Recently the room temperature ferromagnetism 
was reported in the Mn-doped magnetic semiconductors 
(Ga,Mn)N£l which also belong to the III-V family. 

As one of the interesting quantum structures made 
of the DMS materials, the self-organized (In,Mn)As 
quantum dots were successfully fabricated by growing 
(In,Mn)As on the (100), (211)B, (311)B GaAs substrates 
using low-temperature molecular beam epitaxy, which 
were subsequently capped by the GaAs layer .0 The elec- 
tron diffraction pattern and the atomic force microscopy 
measurement confirmed the formation of the (In,Mn)As 
quantum dots. The more detailed experiments after- 
wards revealed that the majority of the Mn atoms were 
actually substituted in the In-site in (In,Mn)As quan- 
tum dots shown by the fluorescence extended x-ca* ab- 
sorption fine structure analysis!! In their samplesJ3'El the 
(In,Mn)As quantum dots were irregularly placed and em- 
bedded in the GaAs. The major challenge in experi- 



ment is to fabricate the uniform and regular array of 
DMS quantum dots and to explore the various physical 
properties of these systems. The interplay between the 
quantum-confined magnetic dots and the non-magnetic 
semiconductors can be very interesting, which has a pos- 
sible application as a tool to imolepient quantum bits, 
a large-scale quantum computer ,EhL3 and other quantum 
devices. 

In this paper, we theoretically study the ferromag- 
netism in diluted magnetic semiconductor quantum dot 
arrays embedded in semiconductors. In contrast to the 
bulk ferromagnetic semiconductors, the following two 
processes lead to the ferromagnetism in the present sys- 
tems: (1) the localized spins of Mn 2+ interact with the 
quantum-confined holes in the dots through exchange 
coupling, and (2) the localized orbital of the quantum- 
confined holes in the dots hybridizes with the itinerant 
holes in the semiconductor valence band. The hybridiza- 
tion allows hole transfer between the DMS quantum dots, 
which may induce the long range ferromagnetic order of 
the localized spins. In turn, the carrier spins both in the 
DMS quantum dots and in the semiconductors are polar- 
ized, which is crucial to realize spintronics. In order to 
describe the basic physics, we propose an Anderson-like 
model Hamiltonian, where the onsite Coulomb interac- 
tion is replaced by the exchange interaction. Our results 
provide a basis for the experiments on exploring the fer- 
romagnetic properties of these systems. 

The paper is organized as follows. In Sec. II, an 
Anderson-type model Hamiltonian is introduced to de- 
scribe the DMS quantum dot arrays embedded in semi- 
conductors. We use the Weiss mean field approximation 
and the self-consistent spin wave approximation, which 
are modified for the present model. In Sec. Ill, we show 
the numerical results of the temperature dependence of 
various physical quantities: the spontaneous magnetiza- 
tion of the localized spins and the spin polarizations of 
the carriers both in the quantum dots and in the semi- 
conductor. Finally we conclude with a brief summary in 
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Sec. IV. 



II. THEORETICAL MODEL AND 
SELF-CONSISTENT APPROXIMATION 

We consider the regular arrays of DMS quantum dots 
with the simple cubic structure embedded in the semicon- 
ductor such as (Ga,Mn)As or (In,Mn)As quantum dots in 
the GaAs Layer. The model Hamiltonian of the system 
can be written by 
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Here Cko- and di a are the fermion operators for the carri- 
ers in the semiconductor and in the DMS quantum dots, 
respectively. For simplicity, we use the parabolic band 
for the carriers in the semiconductor = hk 2 / (2m*) 
with the effective mass to*. Ed is the discrete energy 
level of the carriers within the quantum dots, and fi is 
the chemical potential. The exchange coupling between 
the confined holes and Mn 2+ ion impurity spins is special 
for the DMS quantum dot system and J stands for the 
exchange coupling strength. Sj is the local spin of Mn 
ion impurity and represents the spin of the confined 
holes at the i-th site of quantum dot arrays, which can be 
written as g Yla-a- 1 d\ a f aa r di<j> , where r are the three Pauli 
spin matrices. represents the i-th site of the DMS 
quantum dots in the arrays. The last term in the Hamil- 
tonian (|l|) takes into account the hybridization between 
the holes confined in the quantum dots and the itinerant 
holes in the semiconductor valence band. The present 
Hamiltonian (|l|) is quite similar to that of the periodic 
Anderson model for the heavy fermion compoundstil with 
the onsite screened Coulomb interaction. Here the onsitc 
interaction is replaced by the exchange coupling between 
the localized spin and the hole confined in quantum dots. 

Performing the coarse graining proced-uresu and using 
the Holstein-Primakoff transformations |13 the localized 
spin Si can be written in terms of the bosonic operators 
aj , at as follows 
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where c is the mean number of the magnetic ions Mn 2+ 
in the DMS quantum dots. In the following, the approx- 

a,i ~ V2cS is applied. 
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By transforming from the lattice space to the momen- 
tum space and limiting the involved momentum of the 
carriers in the semiconductor valence band within the 
first Brillouin zone of the corresponding quantum dot ar- 
rays, the Hamiltonian can be rewritten as follows, 

H = E( £k ~ ^) c k<T C k CT + ^2{£da - ^)dl a d ka 
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(!) where e da = e d + f cr with A = JcS, V e $ = V^dP/b 3 , b 
is the lattice constant of the quantum dot arrays, a the 
effective radius of hybridization, and M the total number 
of dot sites. The summations of q, q\, and qi in the last 
two terms in the Hamiltonian (^|) are restricted to the 
values less than the Debye cutoff q c for the spin waves of 
local impurity spins with the relation q 3 , = 6ir 2 c/b 3 . 

In functional integral representation, the partition 
function Z for the Hamiltonian (|^) is given by 
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where the Lagrangian L can be written as 

L = ^2(^1^0^ + dl a d T d ka ) + ^ a q <9 T 
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Here the fermionic and bosonic degrees of freedom are 
represented by the Grassmann variables and the complex 
variables, respectively. 

The itinerant carrier degrees of freedom in the semi- 
conductor valence band can be integrated out easily, and 
we have 



Z = Z, 
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where Z s is the partition function for the free carriers in 
the semiconductor valence band, and L[(fid, d^d] is given 
by 
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where (d T — /z + e ka ) 1 is the Green's function of the 
itinerant carriers in the semiconductor valence band. 
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By subsequently integrating out the remaining carrier 
degrees of freedom in the DMS quantum dots, we finally 
obtain the following partition function 



Z = Z s J £>[a t a]e~' Seff , 



(10) 



where the effective action is given by 
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The mean- field part (G ) is written by 
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and the fluctuation part 8G~ l is given by 
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Hence the exchange coupling and the hybridization of the 
carriers induce the effective non-local interaction between 
the magnetic impurities. 

Expanding Eq. (^) to the quadratic order in and 
a, we obtain 

Sefi = -p ^ a^(w m )£) _1 (q, u/ m )a q (ii/ m ), (14) 
where the inverse of the spin wave propagator is given by 
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and the Green function G^(k, iw n ) is written by 
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The spin wave dispersion is obtained by the following 
analytic continuation iv m — > Q + i0 + , which is given by 
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After the Matsubara frequency summation, the spin wave 
dispersion can be derived as follows 
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The index 



a represents + or — for upper and lower hybridization 
bands, respectively and f(x) is the Fermi-Dirac distri- 
bution function. If the last summation terms in Eq. 
( |l8| ) are ignored, we obtain the Weiss mean field ap- 
proximation with the Weiss mean field given by £1 — 

-J/(2M)Yfaa aA k*f( E ka)> which is proportional to 
the spin polarization of the carriers in the DMS quan- 
tum dots. Here A(q) is introduced as a phenomenologi- 
cal renormalization factor in the last summation terms 
in Eq. (18), which effectively takes into account the 
higher order corrections in the fluctuations in S x and 
SP: A(q) = A + (1 - A ) tanh 2 (g/g c ). Since the Hamil- 
tonian of the system is spin rotationally invariant and 
the ground state possesses a spontaneously brokep_.sym- 
metry, the system has a gapless Goldstone mode.E3 The 
parameter Ao can be determined by imposing the fol- 
lowing condition that f2 q = as q goes to zero. At 
large values of q, the renormalizations in the fluctuation 
term are negligible and so A(q) approaches to one. It 
appeared that for bulk ferromagnetic semiconductors, no 
gap arises in the self-consistent spin wave approximation, 
which we believe is crucial to produce quantitatively re- 
liable resultsfl Hence the modification of the above ap- 
proximation by putting in the renormlization factor A(q) 
can be very important for the present system. 

At finite temperatures, the spin wave dispersion can 
be generalized by imposing the following self-consistency 
conditions for the finite temperature exchange gap A(T) 
given by A(T) = J p d(S z ). Here (S z ) represents the ther- 
mal average of the Mn ion spins in the DMS quantum 
dot arrays, which, are approximately calculated by the 
following formulae! 
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where B s (x) is the Brillouin function and ub (x) the Bose- 
Einstein distribution function. The second term in the 
second equality of Eq. ([H]) describes how the thermally 
induced spin waves from the independent Bosc statistics 
with no limit in the number of spin waves reduce the 
magnetization of the system. The third term takes into 
account the correct spin kinematics, which rules out the 
unphysical states. When fi q is independent of q, for 
exmaple, the Weiss mean field approximation, Eq. ( |2l| ) 
is reduced to the widely used formula (S z ) = cSB s (f3STi). 

On the other hand, using the Hamiltonian (|^) without 
the exchange coupling, one can integrate out the car- 
rier degrees of freedom in the DMS quantum dots to ob- 
tain the Green's function for the semiconductor carriers, 
where the role of the hybridization is taken into account. 
It is given by 



G*(k, iw„ 



-1 



iuJ n - (e k - M) - V es/i ilJ n ~ (^do - M)] ' 

The densities of the carriers in the semiconductors and 
in the DMS quantum dots can be calculated from 
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where V is the volume of the system. The total carrier 
density is given by n tot — ^ a (n* + n^.), which is fixed 
and determines the Fermi energy. 

The spin polarization of the carriers is defined as P s = 
(nj — n|)/(nj +n|) in the semiconductor and Pd — (n^ — 

n f)/( n i m the DMS quantum dots. This definition 

of spin polarizations is different from those commonly 
used in the transport properties, which is defined as_the 
corresponding density of states at the Fermi energyJlj 

Based on the above framework, one can calculate the 
spontaneous magnetization of the localized spins in the 
DMS quantum dot arrays and the spin polarizations of 
the carriers both in the DMS quantum dots and in the 
semiconductor by solving the set of coupled equations in 
Weiss mean field approximation and self-consistent spin 
wave approximations as well. 



III. RESULTS AND DISCUSSIONS 

We have chosen the typical material parameters! of 
the bulk (Ga,Mn)As for the DMS quantum dots, i.e. 
J = 0.15 eV, to* = 0.5to c . The lattice distance b 
of the DMS quantum dots is set to be b — 1.0 nm, 
which is within the spin-coherence length in semiconduc- 
tor GaAs.Ej The carrier bands consist of the upper ej^. 
and lower eZ bands due to the hybridization between 
the discrete energy level of DMS quantum dot and the 
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FIG. 1: Spin polarization of the carriers in DMS quantum 
dots (Pd) and in semiconductors (P 3 ) in Weiss mean field 
theory with the following parameters J = 0.15 eV, m* = 
0.5m e , c* = 0.1 nm -3 , c = 1.0, b = 1.0 nm, V cS = 0.30 eV, 
and e d = 0.46 eV. 




FIG. 2: Spontaneous magnetization of the localized spins in 
DMS quantum dot arrays in Weiss mean field theory with the 
same parameters used in Fig. 1. 



semiconductor valence band. The spin dependence of the 
upper and lower bands is originated from the exchange 
coupling in the DMS quantum dots. Also there exists 
a gap between the spin-up (spin-down) upper and lower 
bands. When the ground state is metallic or semicon- 
ducting, the Fermi level must lie in the spin-down lower 
band. 

In the following we calculate the ferromagnetic prop- 
erties of the present system both in the Weiss mean field 
approximation and in the self-consistent spin wave ap- 
proximation, respectively. 

In Fig. 1 and Fig. 2, we plot the temperature depen- 
dence of the spin polarizations of the carriers in the semi- 
conductor and in the DMS quantum dots, and the spon- 
taneous magnetization for the DMS quantum dot arrays 
in the Weiss mean field theory with the following mate- 
rial parameters c* = 0.1 nm~ 3 , c ~ 1.0, V e g — 0.30 eV, 
and Ed — 0.46 eV. The critical temperature T c for spin 
polarization is about 55 K. One can notice that the spin 
polarizations Pd and P s decrease more stiffly than the 
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FIG. 3: Renormalized spin wave dispersion at 4.2 K in the 
self-consistent spin wave approximation with the same param- 
eters used in Fig. 1. 
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FIG. 4: Spin polarization of the carriers in DMS quantum 
dots (Pd) and in semiconductors (P 3 ) in the self-consistent 
spin wave approximation with the same parameters used in 
Fig. 1. 



FIG. 5: Spontaneous magnetization of the localized spins 
in DMS quantum dot arrays in the self-consistent spin wave 
approximation with the same parameters used in Fig. 1. 
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FIG. 6: Transition temperature T c versus discrete energy 
level Ed in the self-consistent spin wave approximation with 
the following parameters J = 0.15 eV, m* = 0.5m e , c* = 
0.1 ran" 3 , c = 1.0, b = 1.0 nm, and V c ff = 0.30 eV. 



spontaneous magnetization (S z ) of localized spins. The 
T c presumably increases up to the bulk value as one tunes 
the exchange coupling strength J and the hybridization 
V e f[ to higher values. 

Next we go beyond the standard Weiss mean field the- 
ory. In Fig. 3, we have shown the renormalized spin 
wave dispersion curve at T = 4.2 K based on the the 
self-consistent spin wave approximation. The material 
parameters are the same as those used in the Weiss mean 
field calculation. The spin wave frequencies are less than 
the Weiss mean- field value 17 ~ 4.14 meV. As q goes to 
zero, flq oc q 2 , which are correct for ferromagnetic spin 
wave dispersion. As q becomes very large, i7 q approaches 
to the mean field value. Fig. 4 and Fig. 5 illustrate the 
temperature dependence of the spin polarizations of the 
carriers in the semiconductor and in the DMS quantum 
dots, and the spontaneous magnetization for the DMS 
quantum dot arrays in the self-consistent spin wave ap- 
proximation with the same material parameters as be- 
fore. The magnetization and the spin polarizations in 
the self-consistent spin wave approximations drop much 
rapidly comparing to those in the Weiss mean field ap- 



proximation due to the strong spin fluctuations. The T c 
drops to about 22 K less than a half of the mean field 
value of T c , which also happens in the bulk ferromag- 
netic semiconductors. a The qualitative behaviors remain 
the same as the mean field results. 

Fig. 6 shows the dependence of the critical temperature 
T c as a function of the discrete energy level Ed of quan- 
tum dot. The T c increases with the decrease of Ed- As 
Ed decreases, the carriers in the semiconductor easily hop 
to the localized levels and vice versa, which subsequently 
enhances the ferromagnetic coupling between impurity 
spins. The monotonic dependence of T c on Ed is only cor- 
rect in the low carrier density region, in which our model 
is suitable, that is, the Fermi energy level lying in the 
spin-down lower band. In the low carrier density limit, 
the Fermi energy level always lies below e^j_. Obviously 
the enhancement of the hybridization strength between 
the localized carriers in DMS quantum dots and the itin- 
erant carriers in semiconductor makes carrier transfer 
easier leading to the increase of T c . 
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IV. SUMMARY 

We have theoretically studied the origin of the fer- 
romagnetism in diluted magnetic semiconductor quan- 
tum dot arrays embedded in semiconductors based on an 
Anderson-type model Hamiltonian. The hybridization 
between the quantum-confined holes in the DMS quan- 
tum dots and the itinerant holes in the semiconductor 
valence band makes the hole transfer between the DMS 
quantum dots, which induces the long range ferromag- 
netic order of the localized spins in the DMS quantum 
dot arrays through the exchange coupling. Currently 
the available DMS quantum dot systems have irregular 
nanostructures, and the uniform and regular arrays of 



DMS quantum dots are expected in experiments. Our 
results provide a basis for exploring the magnetic prop- 
erties of these systems. 
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